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ABSTRACT 

We present measurements of the FUV (1530A) and NUV (2310A) galaxy luminosity functions (LF) 
at low redshift (z < 0.2) from GALEX observations matched to the 2dF Galaxy Redshift Survey. 
We split our FUV and NUV samples into two XJV—bj color bins and two redshift bins. As observed 
at optical wavelengths, the local LF of the bluest galaxies tend to have steeper faint end slopes 
and fainter characteristic magnitudes M* than the reddest subsamples. We find evidence for color 
dependent evolution at very low redshift in both bands, with bright blue galaxies becoming dominant 
in the highest redshift bin. The evolution of the total LF is consistent with an ~ 0.3 magnitude 
brightening between z ~ and 0.13, in agreement with the first analysis of deeper GALEX fields 
probing adjacent and higher redshifts. 

Subject headings: ultraviolet: galaxies — galaxies: luminosity function, evolution 



1. INTRODUCTION 

The importance of estimating luminosity functions 
(LF) as a function of redshift, color, environment, wave- 
length, etc., to understand galaxy formation and evolu- 
tion, has been well emphasized in the literature. Mul- 
tivariate luminosity functions are crucial to constrain 
theoretical models on small scales where the physics is 
most complex. The shape of the galaxy LFs results from 
non trivial physical processes (Benson et al. 2003, Bin- 
ney 2004) and understanding their evolution is a chal- 
lenging task for numerical simulations (e.g. Nagamine 
et al. 2001). The Galaxy Evolution Explorer (GALEX) 
mission has recently opened a new field of constraints 
by allowing such statistical measures as LFs to be per- 
formed in the rest-frame Ultraviolet (UV) in what is now 
considered the low redshift Universe (z < 1 — 2), where 
most of the evolution of galaxies is thought to have taken 
place. 

In this paper, we use data from the GALEX All-Sky 
Imaging Survey matched to spectroscopic data from the 
2dF Galaxy Redshift Survey to investigate the color 
properties of local UV selected galaxies, and the color 
dependence and evolution of their LF at very low red- 
shift (z < 0.2). Whereas much effort has been devoted 
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to estimate the evolution of galaxy light over the widest 
and furthest range of redshift possible, very few surveys 
have allowed to probe evolution at recent epochs (Love- 
day 2004), and none so far at UV wavelengths. This work 
is a companion paper to Wyder et al. (2004) who present 
the total LFs of FUV and NUV selected galaxies in the 
local (z < 0.1) Universe (hereafter Paper I). It also com- 
plements the work of Arnouts et al. (2004) , Schiminovich 
et al. (2004) and Budavari et al. (2004) who quantify 
the evolution of the UV LF at higher redshift using the 
GALEX Medium and Deep Imaging Surveys. 

The data are summarized in Section 2. Section 3 
presents the color properties of the samples. We es- 
timate and discuss the dependence of the UV LFs on 
galaxy color and redshift at recent epoch in Section 
4. We assumed a flat ACDM cosmology with Ho = 
70 km s^Mpc" 1 , Q M = 0.3 and Q A = 0.7. 

2. DATA 

The GALEX field-of-view has a diameter of 1.2°and 
each pointing is imaged simultaneously in both the FUV 
and NUV bands with effective wavelengths of 1530 A 
and 2310 A, respectively. The GALEX instrum ents 
and mission are des cribed bv iMorissev et all l)2004f) and 
iMartin et all l)2004|) . 

The specific data analyzed in this paper are presented 
in Paper I. In brief, they consist of 133 GALEX All-Sky 
Imaging Survey pointings overlapp ing the 2dF Galaxy 
Redshift Survey l|Colless et al JSOQll) in the South Galac- 
tic Pole region. As the NUV images are substantially 
deeper than the FUV, we used the NUV images for de- 
tection and measured the FUV flux in the same aper- 
ture as for the NUV. The apparent magnit udes were cor- 
rected for foreground extinction using the ISchlegel et ail 
( 1998) reddening maps and assuming the extinction law 
of lCardelli e Faff IpSflh . 

The GALEX catalogs were matched with the 2dF cat- 
alog using a search radius of 6". We restricted the sam- 
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Fig. 1.— Color-color diagrams: FUV - NUV vs FUV - bj for 
the FUV (left panel) and the NUV (right panel) selected samples 
respectively. The thin solid lines are a set of models by Bruzual & 
Chariot (2003) covering a range of star-formation histories, ages, 
metallicities and extinctions. The thick dashed lines are dust free 
models by Poggianti (1997), from Elliptical to Starburst. 



pie to areas with effective exposure times t e ff > 60 sec, 
and removed sources with magnitude errors greater than 
0.4 or contaminated by artifacts from bright stars. We 
also removed overlap regions by restricting the coverage 
to the inner 0.45° of each field. Finally, we excluded 
GALEX sources where the 2dF redshift completeness 
was less than 80% as well as areas of the sky excluded 
from the 2dF input catalog. Our final working area of 
GALEX-2dF overlap is 56.73 deg 2 in both bands. 

We applied a faint magnitude limit of 20 in both 
bands corresponding to the bluest UV — bj color observed 
and beyond which the redshift completeness begins to 
drop. (N.B.: bj was shifted by —0.051 into the AB sys- 
tem. All quoted magnitudes are AB). We also applied a 
bright magnitude limit of 17 to avoid extended sources 
with potential photometric problems. We visually in- 
spected the resulting catalog of 2dF spectra and further 
removed 27 objects with broad emission lines (the ob- 
vious AGN/QSO's). Our final catalogs consist of 1292 
FUV-selected galaxies and 1869 NUV-selected galaxies 
with 17 < UV < 20 and available bj magnitude and 
2dF redshift (with 2dF redshift quality parameter > 3). 
Both redshift distributions extends to z ~ 0.25 (Paper 
I) . FUV fluxes are available for all but one of the galaxies 
in the NUV-selected sample, although uncertainties can 
be quite large for the faintest objects. 

3. COLOR PROPERTIES 

Figure n shows color-color diagrams (FUV — NUV vs 
FUV - bj) for the FUV and NUV selected samples (left 
and right panel, respectively). Typical error bars are 
shown in the upper left corners. The thin lines are a 
set of 42 spectral energy distribution (SED) models by 
Bruzual & Chariot (2003) (hereafter BC), from old el- 
lipticals (in red) to young irregular galaxies, assuming 
a range of star-formation histories (r = 1 Gyr to oo), 
ages (0.1 to 12 Gyr), metallicities (0.2 Z Q to Z@) and 
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Fig. 2. — Model fc-corrections in the bj, FUV and NUV bands. 
The thin solid lines show our set of BC models and the thick dashed 
lines the Poggianti SEDs, as in Fig. HI 



extinctions (Ay — to 1.5). The thick dashed lines are 
dust free SED models by Poggianti (1997). The tracks 
show the color evolution of the various galaxy types from 
z = to 0.25 (roughly counter-clockwise). 

The bulk of both samples is consistent with bona fide 
late-type galaxies, with only a handful of objects dis- 
playing elliptical-like colors. We do not observe any 
source with "{/^-optical excess" colors (bluer than the 
bluest model) such as were observed in the FOCA sam- 
ple (Treyer et al. 1998, Sullivan et al. 2000), suggest- 
ing that the so-called excess was to be found in the es- 
timate of the UV flux of these objects rather than in 
some physical explanation. On the other hand we find a 
small fraction of galaxies with FUV — NUV colors some- 
what bluer than any model (~ 15% in the FUV selected 
sample, ~ 7% in the NUV selected sample). This is 
partly due to newly measured, small but opposite offsets 
in the FUV and NUV calibrations, adding up to +0.13 in 
FUV-NUV (new in-flight calibrations, P. Morissey, pri- 
vate communication), not applied in the current catalogs. 
But shortcomings in the models at UV wavelengths can- 
not be ruled out. Despite the wide coverage in parameter 
space, the data are not fully accounted for (e.g. the red- 
dest sources in both FUV-NUV and FUV~b 3 ). Fig.[T] 
suggests that UV selected galaxies may require alterna- 
tive model ingredients (e.g. more erratic star-formation 
histories as suggested by Sullivan et al. 2003). 

We assign a best fit SED (drawn from the combined 
and interpolated BC and Poggianti libraries) to each 
galaxy using its redshift and two independent colors 
(FUV - NUV and FUV - bj). Figure H shows the 
model fc-corrections in the bj , FUV and NUV bands (thin 
solid lines for BC and thick dashed lines for Poggianti). 
UV fc-corrections are small enough for the late-type, low- 
rcdshift galaxies dominating the sample (as opposed to 
fc-corrections in optical bands) that uncertainties on the 
assigned SEDs have negligible impact on the resulting 
LFs. 

We can also fit a power-law through the 2 band- 
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TABLE 1 

SCHECHTER FUNCTION PARAMETERS (0* IS IN MPC -3 ). 



Band 


z 


Type 


M* 




—a 




log</>, 




FUV 


0.0-0.1 


All 


— 18.04 ± 


.11 


1.22 ± 


.07 


—2.37 ± 


.06 




0.1-0.2 




-18.31 ± 


.06 






-2.33 ± 


.07 




0.0-0.1 


Blue 


-17.89 ± 


.15 


1.29 ± 


.09 


-2.54 ± 


.09 




0.1-0.2 




-18.44 ± 


.09 






-2.73 ± 


.10 




0.0-0.1 


Red 


-18.19 ± 


.22 


1.10 ± 


.17 


-2.80 ± 


.12 




0.1-0.2 




-18.15 ± 


.08 






-2.50 ± 


.10 


NUV 


0.0-0.1 


All 


-18.23 ± 


.11 


1.16 ± 


.07 


-2.26 ± 


.06 




0.1-0.2 




-18.53 ± 


.05 






-2.27 ± 


.05 




0.0-0.1 


Blue 


-18.04 ± 


.16 


1.25 ± 


.10 


-2.48 ± 


.09 




0.1-0.2 




-18.71 ± 


.08 






-2.82 ± 


.08 




0.0-0.1 


Red 


-18.28 ± 


.IS 


0.97 ± 


.11 


-2.56 ± 


.09 




0.1-0.2 




-18.36 ± 


.06 






-2.37 ± 


.06 



passes and derive a spectral index for each galaxy (j3 = 
(FUV-NUV)/0A25-2, assuming F\ oc X 13 ). The mean 
spectral indices of the FUV and NUV selected samples 
are (/?) ~ —1.6 and —1.3 respectively. The UV slope 
f3 has been shown to be a good tracer of dust attenua- 
tion in galaxies undergoing a strong starburst, but not so 
reliably in others (Bell 2002, Kong et al. 2004). In par- 
ticular, Buat et al. (2004) found that local NUV selected 
galaxies tend to have lower attenuations than starburts 
for a given slope. As our two samples exhibit slopes cor- 
responding to low attenuations in pure starburst galax- 
ies, we conclude that they must be very little affected by 
dust. (See Paper I for a quantitative discussion). 

4. COLOR AND Z DEPENDENT LUMINOSITY FUNCTIONS 

We split both samples into two rest-frame color bins: 
(FUV — bj)o < and > 2 (roughly corresponding to an Sb 
spectral type and to the mean of the color distribution) , 
referred to as blue and red subsamples respectively. We 
further split each of these subsamples into two redshift 
bins: z < 0.1 and 0.1 < z < 0.2, referred to as low 
and 'higher' redshift subsamples respectively. The mean 
redshifts in these two bins are 0.055 and 0.12 respectively 
for the FUV sample, and 0.057 and 0.13 for the NUV 
sample. 

As in Paper I, the redshift completeness is defined as 
the ratio of UV objects matched to a 2dFGRS counter- 
part with (high quality) redshift to the total number of 
UV galaxies in a given magnitude bin, as computed by 
Xu et al. (2004) using the SDSS star/galaxy separation 
criteria in overlap regions. It is roughly constant over the 
range 17 < UV < 20 in both bands, with mean values 
of ~ 79% in the NUV and 92% in the FUV (cf Fig. 1 in 
Paper I). We computed binned lumi nosity functi ons us- 
ing the traditional Y max estimator l|Feltenlll976j) in the 
12 subsamples (blue, red, and total at low and higher z 
in both bands), and derived best fit Schechter functions 
l)Schechterlll976j) for each of them. Since the magnitude 
range in the highest redshift slice is too bright and nar- 
row to constrain the slope meaningfully, we fixed it to 
the low redshift value in each case. 

Results are shown in Fig. [3] and the best-fit Schechter 
function parameters are listed in Table ^ The FOCA 
2000A LF (Sullivan et al. 2000) at a mean redshift 
z ~ 0.15 is also plotted for comparison (converted to the 
AB magnitude system and to our Hq value but without 



accounting for the difference in bandpasses). We refer to 
Paper I for a discussion on the discrepancy between the 
FOCA and GALEX LFs. As observed at optical wave- 
lengths (Madgwick et al. 2002, Nakamura et al. 2003), 
the local LFs of the bluest galaxies tend to have steeper 
faint end slopes and fainter characteristic magnitudes M* 
than the reddest subsamples (although our definition of 
red would qualify as blue in optically selected samples). 
The fact that the latest types are observed to fainter ab- 
solute magnitudes than the earlier types induces a small 
bias in the estimate of the faint end slope of the full sam- 
ple. The dotted lines in the upper panels of Fig. 01 show 
that the sum of the blue and red LFs is slightly steeper 
at the faint end than estimated for the galaxy popula- 
tion as a whole, due to the shortage of red galaxies in 
the faintest bins. 

The total FUV and NUV LFs in the highest redshift 
bin (lower panels in Fig-EJ are both consistent with ~ 0.3 
magnitude brightening in M* with respect to the LFs at 
z < 0.1. This evolution is quite similar to that found in 
optically selected samples at the same redshifts (Loveday 
2004). Evolution can also be inferred from the steeply 
rising V/V ma x distribution of both z-unlimited samples 
(with (V/V m ax) = 0.54 in both bands instead of the 
0.5 expected from a non evolving population). As the 
redshift distributions show two strong peaks around 0.06 
and 0.11 (Paper I), the appearance of evolution could be 
attributed to large scale structure effects, but removing 
these two structures from the samples did not affect our 
results. 

Evolution also seems to occur predominantly in the 
blue subsamples, with blue galaxies dominating the 
bright end of the LF in the highest z bin. We checked 
that no color trend was observed as a function of bj, as 
might be expected if the 2dFGRS was strongly biased 
against, e.g., faint blue low surface brightness galaxies 
(LSBG), making bright blue galaxies seem dominant at 
higher redshift. The issue of missing galaxies in the 2dF- 
GRS, in particular LSBGs, has been addressed in a num- 
ber of papers, most recently by Cross et al. (2004). The 
authors find that the survey does indeed miss ~ 15% of 
the galaxy population, of which 6% only are identified 
as LSBGs. If this particular class of objects is unlikely 
to affect our results, the remaining missing galaxies may 
bias our analysis, should their redshift distribution look 
different from the one actually measured. A truly UV- 
selected spectroscopic survey is underway to address this 
question. 

The present results are consistent with the FUV analy- 
sis of a much deeper GALEX field probing adjacent and 
higher redshifts (Arnouts et al. 2004, Schiminovich et 
al. 2004). They are also in rough agreement with the rate 
of evolution derived from a photomctric-redshift study of 
GALEX/SDSS data at low z (Budavari et al. 2004), al- 
though our FUV LF in the range 0.1 < z < 0.2 seems 
brighter than the photo-z estimate at overlapping z. 
Both analysis agree very well in the NUV however (we 
postpone attempts at explanations for the FUV discrep- 
ancy until more comparable datasets are available). We 
note that evolution seems to be limited by the observed 
galaxy number counts, especially in the FUV band (Xu 
et al. 2004), so that the rate of evolution detected at 
z < 1 in the GALEX data is expected to slow down. 
This appears to be the case based on a photometric- 
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Fig. 3. — Top panels : The FUV and NUV LFs at z < 0.1 (left and right panel, respectively). The blue and red lines (and associated 
filled circles and asteriks) represent the blue and red subsamples respectively. The black solid lines (and associated empty circles) show the 
LFs of the full samples (Paper I), while the dotted lines are the sum of the blue and red Schechter functions. The green dotted line is the 
FOCA 2000A LF at z ~ 0.15 (Sullivan et al. 2000). Bottom panels : The FUV and NUV LFs at 0.1 < z < 0.2. Same symbols and lines 
as above. 



redshift analysis of the HDF North and South probing 
1.75 < z < 3.4 (Arnouts et al. 2004). 
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